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Abstract 



A physical analysis and mathematical formulation of the structure of an accretion shock 
wave in the case of young stars are given. Some results are presented, and the dependence 
of the structure of the shock on various parameters is investigated. In the general case, 
the relative intensities and profiles of lines in the spectra of T Tauri stars should depend 
on the velocity and density of the infalling gas, as well as on the accretion geometry. The 
calculation results may serve as a base for elucidating the nature of the activity of young 
stars. 



1 Introduction 

T Tauri stars are young (t < 10 7 years), low- mass (M < 2M Q ) stars at the stage of 
gravitational contraction toward the main sequence. From the observational point of view, 
they are late-type subgiants with emission lines of hydrogen and Ca II, as well as of some 
other atoms and ions, mostly with relatively low ionization potentials. T Tauri stars are 
also characterized by excess continuum emission at short visible wavelengths and in the 
infrared compared to main-sequence stars of the same spectral types (see the reviews 
[1, 2]). The variety and relative intensities of their emission-line spectra, as well as their 
continuum excesses, correlate with the equivalent width of the H a line; thus, we may 
consider this line to be an indicator of the activity of low-mass young stars. It is accepted 
to call stars with Wn a > 5 — 10 A classical T Tauri stars (CTTSs), and precisely such 
stars will be the subject of our studies below. 

Infrared and radio investigations of CTTSs testify that these stars are surrounded 
by disklike gas-dust envelopes. IUE observations of CTTSs show the presence of lines of 
highly ionized atoms, up to N V; the luminosity in these lines is several orders of magnitude 
higher than for the corresponding solar lines [3]. However, at X-ray energies, between 0.5 
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and 4.5 keV, CTTSs have luminosities < 10 31 erg/s, i.e., 2-3 orders of magnitude lower 
than expected by analogy with the Sun based on the observed intensities of, e.g., the C IV 
A1550 line. The absence of forbidden FeX A6375 and FeXIV A5303 coronal lines in the 
spectra of CTTSs suggests that this is not associated with the absorption of X-rays in the 
extended circumstellar envelopes, but rather is due to a real difference of the observed 
Lx/Lciv ratio from the solar ratio. The comparatively hard spectrum (T > 3 x 10 6 K) 
and variability of the observed X-ray emission suggest that it is associated with flare 
activity [2, 4]. 

There is not yet any consensus about the nature of the regions in which the CTTS 
emission-line spectrum forms, nor about the origin of the rather strong outflows of matter 
from the vicinitiy of CTTSs in the form of so-called CO flows and jets. Up until the early 
1980s, it was thought that the activity of CTTSs was the result of powerful chromospheres 
and/or coronae, to some extent similar to the corresponding solar structures, and that the 
reason for the existence of hot layers above the photosphere and of the stellar wind was the 
dissipation of mechanical energy generated by the stellar convective zone [5-9]. Later, it 
was proposed that the emission-line spectrum was associated with heating of material from 
the disklike circumstellar envelope falling onto the young star. If the inner edge of the disk 
is in contact with the stellar surface, the heating should take place in a boundary layer, 
where the disk material is decelerated from the Keplerian velocity to the stellar rotation 
velocity [10, 11]. However, if the young star has a fairly powerful dipolar magnetic field, 
then, by analogy with X-ray pulsars, it may be that the inner disk boundary does not 
reach the stellar surface, and that the disk material one way or another penetrates the 
magnetosphere and slides down along lines of force. In this case, it is natural to associate 
the region of formation of the emission-line spectrum with the shock wave that should 
arise during the collision of the infalling gas with the stellar surface [12-14]. In accretion 
models, the observed outflow of matter is interpreted in terms of a magnetohydrodynamic 
wind from the disk surface or directly from the stellar magnetosphere [2, 12, 15]. 

In principle, investigations of the profile shapes and relative intensities of CTTS 
emission lines could elucidate the nature of the hot regions in young stars. However, 
numerous attempts to use lines of atoms and ions with low ionization potentials (e.g., the 
hydrogen Balmer lines) in such studies have not yet been successful: similar profiles are 
obtained for quite different assumptions about the geometry, velocity field, and physical 
conditions in the emission-line formation region [1,2,16,17]. Moreover, high-angular-resolution 
spectral observations have shown that the profiles of these lines are formed both near the 
stellar surface and in the extended circumstellar envelope [18, 19]. 

Analyses of the line emission of ions with ionization potentials above 20 eV are much 
more promising in this connection [20]. At the present, the following is known about the 
profiles of these lines. 

(1) The Hell A4686 lines in CTTS spectra have equivalent widths of less than 1 A, 
and it became possible to investigate their profiles with sufficiently high signal-to-noise 
ratio only recently [21-23]. In all cases (nine stars), the profile of this line is asymmetric 
and redshifted by a few dozen kms -1 . 

(2) In [24], Gomez de Castro et al. analyze the profiles of the resonant lines of the 
Si III, Si IV, C III, and C IV ions using archival high-resolution IUE spectra of four CTTSs 
(T Tau, RU Lup, DR Tau, TW Hya) obtained in the 1980s. The profiles of lines that could 
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be reliably distinguished against the noise background were asymmetric and redshifted. In 
the case of the C IV A1550 doublet lines, two peaks of comparable intensities were reliably 
identified, with their maxima shifted by approximately +50 and +250 kms -1 with respect 
to the laboratory wavelength. 

(3) In the spectra of many CTTSs, there are He I emission lines that are about an 
order of magnitude more intense than the Hell A4686 line. According to [1], the He I 
line profiles have a broad (FWHM ~ 150 kms" 1 ) peak that is symmetric with respect 
to the laboratory wavelength. In a number of cases, a narrower (FWHM < 50 kms" 1 ) 
central component is superimposed. On the other hand, in high-resolution spectrograms 
with good signal-to-noise ratio, some stars have a broad absorption feature in the red 
wing of the A5876 triplet; in RY Tau and BM And, the entire A5876 line is in absorption, 
not emission [25]. 

In [24, 25] , it was proposed that the above features of the lines of helium and ions with 
ionization potentials above 20 eV can be understood if these lines form in a shock wave 
near the surface of the young star. Note that the observations currently available do not 
exclude the possibility that CTTSs have a comparatively powerful chromosphere and/or 
corona, which may contribute to the line emission of ions with high degrees of ionization. 
Therefore, studies of the emission of an accretion shock wave should also elucidate whether 
or not powerful chromospheres are present in low-mass young stars. 

After calculating the shock-wave structure, we can find the intensities of specially 
selected "signal" lines; comparisons with observations enable us to determine the velocity 
and density of the accreted gas, and also the total accretion rate and the area of the stellar 
surface on which the material falls. To determine the "signal" line profiles, it is necessary 
to appropriately sum the emission from different points in the shock wave, so that analysis 
of the observed profiles and, in particular, their variation with the axial stellar rotation 
should yield information about the geometry of the accretion region. 

The structure of radiative shocks has been investigated in a number of studies in 
connection with various problems in the physics of the interstellar medium [26, 27] and 
processes in the envelopes of nonstationary stars [28]. Therefore, we should explain the 
peculiarities of the study of such shocks in the case of CTTSs, and why we cannot apply 
already available results of gas-dynamical calculations to the interpretation of the spectra 
of young stars. The difference from the case of the interstellar medium is due to the fact 
that we consider a shock in a rather dense gas in which the velocity of the stream flowing 
into the front is 150 — 450 kms -1 . The much higher density means, for instance, that the 
forbidden lines of 01, Oil, Nil, etc., cannot play an important role in cooling the gas. 
Moreover, in an interstellar medium with velocities higher than 100 kms" 1 , postshock 
radiative cooling does not play a significant role, and the gas relaxes adiabatically [26]. 
On the other hand, shock waves in the envelopes of nonstationary stars are characterized 
by higher densities and velocities for the stream flowing into the front that are a factor 
of 2-3 lower [28]. Thus, the situation of interest to us is, in some sense, intermediate, and 
has not previously been investigated to our knowledge. 

Before addressing the problem of the structure of an accretion shock wave in the 
case of a low-mass young star, we made some physical assumptions reflecting the specific 
character of the situation investigated, which made it possible to considerably simplify 
the general system of radiative gas-dynamical equations. Therefore, the main goal of this 
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work is to show that the results of our numerical calculations do not contradict our initial 
assumptions, and can serve as a basis for the quantitative interpretation of the spectra of 
young stars. In the last section, we will present data on the structure of accretion shocks 
and the calculated intensities of some observable resonant lines. 

2 Formulation of the problem 

Consider a young star with mass M* and radius with a large-scale dipolar magnetic 
field surrounded by an accretion disk. Suppose that the stellar magnetic field stops the 
accretion disk at some distance Ri > R*; the value of Ri depends on the accretion rate 
M a and the intensity of the stellar magnetic field [12-14]. The disk material is frozen in 
the magnetic lines of force and slides along them under the action of the force of gravity, 
acquiring near the stellar surface a radial velocity [29] 

(Here, we assume that the angular velocity of the stellar rotation is much lower than the 
Keplerian velocity (uk = (GM^y^R* 3 ^ 2 ). According to the data of [30] and [25], for the 
stars RU Lup and RY Tau, V ~ 170 and 400 kms -1 , respectively. Therefore, we will be 
interested in velocities for the infalling gas from 150 to 450 kms -1 . 

During the collision of the gas with the stellar surface, a shock wave should arise. Let us 
assume that this shock has the following qualitative structure [30]. The gas compressed in 
the shock front is strongly heated, then gradually decelerates and cools to a temperature of 
the order of the stellar effective temperature, smoothly passing into the stellar atmosphere, 
which is in hydrostatic equilibrium. The cooling takes place via volume losses to emission 
in resonant lines of the ions of the most abundant elements; the largest fraction of the 
energy is in photons with energies above 13.6 eV. Half of the photons from the cooling 
zone move away from the star and are absorbed by the preshock matter, creating a heating 
and ionization zone that is transparent in the continuum at A > 912 A. The other half of 
the photons are absorbed in the transition region between the radiative-cooling zone and 
the stellar atmosphere, and are then reemitted in the continuum. Only a small fraction 
of these photons goes into the Lyman continuum, so that the emission of the transition 
region and, consequently, its structure, virtually do not influence the structure of the 
higher-lying regions. 

Within the framework of these assumptions, we estimated in [30] the number density 
N (cm -3 ) of the accreted gas far upstream from the front: it turned out that situations 
in which 10 < logiVo < 13 are of the most interest. For this interval of densities, we will 
also assume that the typical dimensions of the postshock radiative-cooling region and of 
the preshock HII region are smaller than the stellar radius, and are also smaller than 
the extent of the front perpendicular to the gas motion. Therefore, we will consider the 
shock wave to be plane-parallel and will solve the corresponding one-dimensional problem. 
Finally, we assume that the shock front is at rest with respect to the stellar surface, and 
that all the shock parameters depend only on the spatial coordinate. Thus, we consider 
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below a solution of the one-dimensional, radiative, gas-dynamical equations describing 
the structure of an accretion shock wave that is stationary in Eulerian coordinates. 

3 System of radiative hydrodynamical equations 

In a stationary, one-dimensional shock wave, all parameters depend only on the spatial 
coordinate z, and it is assumed that (d/di) z = everywhere in the gas-dynamical 
equations. We will also assume that a chaotic magnetic field is frozen in the infalling 
plasma, and that this field is much weaker than the large-scale stellar field channeling the 
flow of matter, but is still sufficiently intense to suppress the electron thermal conductivity. 
We are interested only in the field component perpendicular to the z axis; we designate 
the intensity of this component H. 

We will take the accreted matter to consist of a mixture of H, He, C, N, O, Ne, Mg, 
Si, S, and Fe, with the relative abundances 

N a ^ N a 

a 

where N a and iV are the volume number densities of a given element [31] and the total 
number density of the nuclei of all elements, respectively. We denote Nf and x\ to be the 
volume and relative number densities of ions of the ith degree of ionization of element a 
(0 < % < Z a , where Z a is the nuclear charge of this element), so that 

^ a = f>-, f>? = l. 

i=0 i=0 

Let N e be the volume number density of free electrons. Then, designating x e = N e /N, we 
can write the equation of electric neutrality of the plasma 

x e = xf + xf c + 2x% c + 3 x 1(T 4 , (1) 

where the last term takes into account the single ionization of C I by lines of the hydrogen 
Lyman series in preshock HI regions. Note that, at this stage, we are interested solely in 
lines of ions with high ionization potentials (see above). Therefore, we will not investigate 
the structure of the H I regions in detail; we will show later, however, that these regions 
virtually do not affect the fluxes in the lines of interest to us. 

In the absence of molecular hydrogen, the continuity equation is [32] 

j = NV = const = N V , (2) 

Here and below, all quantities with the subscript refer to the preshock unperturbed 
flow, i.e., formally to z — — oo. We will assume that the gas temperatures for the various 
ions and atoms are the same and are equal to Tj, but, generally speaking, are different 
from the temperature T e of the electron gas. We will express the gas density in terms of 
the mean molecular weight ^ : 

p = J2 ™ a N a = fi imu N, & = CaA a ~ 1.3, (3) 
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where m a and A a are the mass and atomic weights of an element. Since the pressure of 
an ideal gas is P = NkT, the momentum conservation law in our case can be written 



H 2 H 2 
NkTi + N e kT e + famuNV 2 + — = const ~ N kT i0 + ^mnN V 2 + (4) 

assuming that the gas is weakly ionized at z — — oo. Another algebraic equation describing 
the flux freezing of the small-scale chaotic magnetic field is [33]: 

- = — (5) 

We will now consider the energy- conservation equations. In the two-temperature approximation 
adopted, separate equations must be written for the electrons and heavy particles, i.e., 
atoms and ions. For the heavy particles, the energy sink is the exchange of energy with 
electrons, and therefore the energy conservation law for them per heavy particle is [32] 

v i (H + VNm i ih) = -"" N - (6a) 

Here, we have used the stationarity condition, which allows us to replace ^ with V (J^) . 
Energy exchange with ions is much more efficient than with neutral atoms. Since, for 
V ~ 300 kms -1 , the most abundant element - hydrogen - is appreciably ionized over 
nearly the entire flow region of interest to us, we will take into account in (6a) only 
interactions between electrons and protons. In this case [32], 

3 jj Ti — T e \k 

u ei ~ -k^x, -^pr^ 

where \ K ~ 9.4 + 1.51nT e — 0.5 In N e is the so-called Coulomb logarithm. 
For free electrons, the energy conservation law per heavy particle is [32] : 

{^-x e kT)j + VN ekT e ^ (J^j = Q ph - Q col N e + u ei N e . (7) 

The quantity Q p h describes the heating of gas by photoionization; Q co \ denotes the sum 
of terms associated with energy losses to excitation Q exc and ionization Q ion via electron 
collisions, and with radiative cooling during recombinations Q rec . In other words, 

Qcol Qexc + Q ton + Q rec- 

We will now write these terms in explicit form. 

We assume that the cooling postshock gas is transparent to its own "cooling" radiation 
down to a temperature of ~ 10 4 K, and that the heating and ionization of this gas by this 
radiation takes place upstream from the shock front and far downstream from the front, in 
the so-called transition region. We can write for the heating of gas by "external" radiation 

Gi* = E^-ECEs?E vUhn - xy , (8) 

l nU l a i i 
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where Fi, is the flux of radiation in the spectral line with frequency v\\ xij i s the 
ionization potential of the i-th ion of element a from shell j; afj t is the cross-section for 
photoionization of this ion by the line I; and the numerical value of v\\ Xij was calculated 
using formulas from [34]. (We also described the ionizing continuum radiation using certain 
fictitious lines - see below). We considered photoionization from the inner K and L shells 
under the assumption that this is always accompanied by the Auger effect, excluding H-, 
He-, and Li-like ions. In equation (8), we took this into account by substituting Ixth m 
place of XiK- Since the flux of photons with energies above 1 keV is small for the velocities 
V of interest to us, we neglected the possibility of ionization from M shells. 

For a normal abundance of heavy elements, energy losses due to collisional ionization 
are significant only for hydrogen and helium, so that 

D — t ^ H /r H -v H _l_ t ^ He /r He -v He _l_ t ™He He -He / Q \ 
Qion - CH^o % Xo + CHe% % Xo + CHc^l <?l Xl , (»J 

where qf is the coefficient for ionization by electron collisions. Energy losses via recombination 
were also taken into account only for hydrogen and helium. If an ensemble of free electrons 
loses an energy /3(T)kT to recombination, then (cf. [35]) 

Q rec = kT (Ch^Xo + (Hex? e f3? e a?° + ( He x? e f3* e a?°) , (10) 

where i is the coefficient for radiative recombination for an ion with charge % + 1 to 
become an ion with charge %. We assumed that the optical depth in the resonant lines 
of hydrogen and helium is fairly high, so that only recombinations to excited levels were 
taken into account in (10). 

We will now present an expression for the energy losses to excitation of discrete 
levels, with the subsequent emission of the corresponding photons. We considered electron 
collisional excitation only from the lowest term of the ground configuration, so that 

Q exc = E Ca E < E h ^ (*iCZS - cm. (ii) 

a i k 

where n k is the relative population of level k, C^' l k and C k \ are the coefficients of electron- 
collisional excitation and deexcitation, respectively, and e\ t k is the probability of escape 
of a photon from the emission region, which we assume to be equal to 1. We took our list 
of the most intense lines (in total about 150) from [36], and the expressions for C"'^ from 
the reviews [37]. 

Consider the ionization equilibrium equations. Let J 4 a be the rate of photoionization of 
the ith ion of element a from all its subshells per unit volume. Then, assuming afj t = 0, 
if hv\ < Xij, we obtain 

/? = ^^^^ 

i 1 j 

for i < Z a and If = for i = Z a ; the subscript j denotes the electron shells of the 
ion considered (1 < j < j max )- J? is an analogous quantity characterizing the rate of 
photoionization only from those shells for which the removal of an electron is accompanied 
by the Auger effect. 

Generally speaking, the transition of an ion % to the i — 1 state may take place as a result 
of radiative, dielectronic, or three-particle recombination. In the situation considered 
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here, nearly everywhere logiV e > 10. According to [38], at such densities, collisions 
with electrons disrupt highly excited ionic states, leading to a substantial decrease of 
the dielectronic recombination rate, since captures to high-lying levels play the major 
role in this process. For logiV e < 15, the radiative recombination rates depend much 
more weakly on the density, since the contribution of high-lying states is much smaller 
in this case. Using the data in [39-41], we came to the following conclusions about the 
dielectronic recombination rate for the interval of densities and temperatures of interest 
to us: (1) for ions with charges smaller than six, it is considerably lower than the radiative 
recombination rate; (2) for hydrogen-and helium-like ions, density-dependent corrections 
do not exceed 20%. Therefore, for the CII-CIV, NI-NV, OI-OVI, Nel-NeVI, Mgll- 
MgVI, Si II Si IV, SII-SVI, and Fell-FeVI ions, we set the dielectronic recombination 
coefficient to zero, while, for the remaining ions, we neglected the density dependence, 
as for the radiative recombination coefficient. In Section 5, we will present arguments 
indicating that this approximation is unlikely to introduce significant errors into the final 
result. 

Three-particle recombination plays an appreciable role only at rather high densities; 
the main contribution to this process comes from captures to high-lying levels. In particular, 
for hydrogen at a temperature of 10 4 K and with log N e = 14, the rate of three-particle 
recombination to the first level is approximately three orders of magnitude smaller than 
the total radiative recombination rate [32]. Interactions with ions and electrons in a plasma 
with such parameters limit the number of stable levels of the hydrogen atom to a value 
not greater than 100 [42]. Apparently, this means that, in the problem at hand, the 
rate of three-particle recombination for hydrogen nowhere exceeds ~ 10% of the radiative 
recombination rate; for this reason, we did not take it into account. Since the rate of three- 
particle recombination rapidly decreases with increasing nuclear charge [43], we likewise 
did not take it into account for the remaining ions. 

We denote the radiative and dielectronic recombination coefficients af i _ l and 
respectively. The number of recombinations involving the ith ion of element a in unit 
volume per second can be written 

it- _ x = N e + 

for i = 1 — Z a and 

Rh-i = 

for % — 0. 

In terms of this notation, the ionization-equilibrium equations take the form (for 
simplicity we omit the subscript a): 

dx ; 

V-^- = Ji- 2 Xi-2 + (qi-iN e + - Jj_i) - (qiN e + k + R iti -i) + R i+lti x i+1 . (12) 

We calculated a(T), k(T), and q(T) using formulas from [36]. For a number of ions, charge 
exchange reactions of the form 

A +n + H _^ A +(n-l) + H + 

A +n + He ->■ A +in - 1] + He + 
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play an important role in the ionization balance. Formally, we took this into account by 
adding the terms N^x^x^b^iT) and/or iV CiieXiX^b^'^ \T) to equation (12) if the rate 
of the corresponding reaction (6"' H or 6"' Hc ) exceeded 10~ 14 cm -3 . We also considered the 
reverse reactions with H + in accordance with [44] for the atoms N° and 0°. 

We will assume that ionizing radiation is generated only behind the shock front, and 
that the generation region has a small optical depth in all lines. In this case, we can write 
for the flux of radiation propagating away from (or toward) the star 

The expression for the flux in lines with energies below 13.6 eV forming behind the front, 
but outside the transition region, has a similar appearance. Below, we will ascribe the 
subscript I to lines with energies above 13.6 eV and use the subscript k in the general 
case. 

In [30], we showed that, in the problem at hand, the temperature of the postshock 
electron gas is appreciably lower than 10 7 K. At such temperatures, the gas is cooled 
primarily by discrete transitions, but the continuum radiation makes a significant contribution 
to the X-ray emission [36]. Therefore, when calculating the spectrum of the cooling zone, 
we took into account recombination emission by H, He, and the C VI, C VII, O VIII, and 
OIX ions [36], as well as free-free emission [45]. In contrast to [36], we did not take into 
account two-photon emission, since its contribution is small when logiV e > 11 [46]. 

If self- ionizing radiation is not generated upstream from the front, the variation of the 
flux of ionizing photons with distance from the front (z = 0) is described by the law 

F,(z) = F,(0)exp(-7i), (14) 

with 

where the summation is over elements, ions, and shells for which hv\ > xtj- Jumping 
ahead for a moment, we note that the contribution of lines with A < 912 A to the cooling 
function of the preshock gas proves to be insignificant, which justifies a posteriori the use 
of approximation (14). 

The kinetic energy of the directed motion of the gas is converted into the thermal 
energy of atomic nuclei in a layer that we call the shock front. According to [32], its 
thickness is a factor of (mp/mg) 1 / 2 ~ 40 smaller than that of the region following it, 
in which the temperatures of the electrons and heavy particles are equalized. We will 
treat the shock front as a mathematical discontinuity at z — 0, and ascribe the subscript 
1 to quantities immediately upstream from the jump and the subscript 2 to quantities 
immediately downstream from the jump. For heavy particles, the relationship between the 
parameters upstream and downstream from the jump is obtained by solving the system 
consisting of equations (5)-(8) and the energy conservation law: 

ArT , (nimnV? 5kT a \ V X H\ Ar T , //i*m H V; 2 5kT i2 \ V 2 H\ 
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Introducing the Alfven velocity V a = H/(ATTfi i muN) 1 ^ 2 and designating 



V 2 
v Al 



Vi 



2 ' 



B 



2 ' 



C=l + ^ + 55, 



we obtain 



N 2 



Va2 = Vai\ tt, 

V Vl 



8A 



rp _ Til 



Vo 



A 



V 2 AV 1 



V 2 2V 2 



(16a) 
(166) 
(16c) 



In the case of a strong shock and a weak magnetic field, i.e., when A, B — > 0, we obtain 
from (16) the well-known relationships: 



N 2 



Yi 

Vx 



1 

I' 



i2 



(16d) 



Recall that, if we use the one-temperature approximation everywhere, then, in the strong- 
shock limit, the temperature immediately behind the front will be a factor of two lower 
than that indicated by formula (16d). The electron gas at the viscous jump is compressed 
adiabatically, and the degree of ionization remains virtually unchanged [32]. Hence, we 
obtain 

,y.2/3 



Vo 



(16e) 



In the strong-shock limit, we have T e2 ~ 2.5T el . 



4 Boundary conditions and the solution method 

Before proceeding to a direct analysis for the various flow zones, note that, by adding (6a) 
and (7), we can obtain with (4) the relationship 

\k (Tj + x e T e ) + m H /^ + 4^ " Ch*K - Chc^ c (x? c + X?j - Chc*?^? 

= -2 ^ F k - 2 £ F- c (hu < x h) , (17) 

k+l H,He 

expressing the law of conservation of energy for the problem at hand. The term in square 
brackets is the difference of the energies of the infalling gas far upstream and downstream 
from the shock, while the right-hand side contains the energy flux carried away by photons 
with E < 13.6 eV. Photons with E > 13.6 eV are absorbed upstream from the shock front 
and in the upper layers of the stellar atmosphere, and are then reprocessed into radiation 
with E < 13.6 eV, which escapes from the system. 

The gas temperature T far upstream from the front should be of the order of several 
thousand Kelvin, i.e., ~ 0.5 eV. The magnetic energy associated with the infalling plasma 
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is even smaller (see below). On the other hand, for Vq = 300 kms -1 , E^m = /-t^H^ 2 /2 — 
550 eV. It therefore follows from (17) that the energy of the infalling material is nearly 
entirely in the form of kinetic energy, which is spent primarily in radiative cooling, while 
the losses to gas ionization are small. It is convenient to express the line emission flux 
-Ffc(O) in units of jEki n /2, introducing the quantities 5k characterizing the fraction of 
energy emitted in the line k downstream from the shock front: 

Ft(0) = (18) 

with J2kh < 1. 

4.1 The preshock region 

In accordance with [26, 32], we expect that the velocity and density of the preshock gas 
change little, so that 

V = V {1-5V), N = N (1 + 5N), H = H (1 + 5H), (19) 

where 5V, 5N, 5H < 1; as a consequence of (2), 5V = 5N, and it follows from (5) that 
5H = 5N. 

— i — i — i — | — i — i — i — i — | — i — i — i — i — | — 




J i i i i I i i i i L_ 

5 10 15 

T 



Phc. 1: Electron (T e ) and ion (Tj) temperalures as functions of the optical depth r in 
a hydrogen plasma (see text). The thick line shows the run of the temperature when 
T e = T t = T. 

Upstream from the shock front, the electrons are heated directly by photoionization, 
and we cannot say in advance whether the energy exchange with heavy particles will 
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be efficient enough for their temperature to become essentially the same as that of the 
electron gas. We can answer this question by solving a model problem in which pure 
hydrogen is heated by a single ionizing line and the preshock gas has a constant velocity 
V - V . 

We will provisionally neglect recombination, ionization, and radiative losses due to 
electron-collisional excitation. We then obtain from equations (6)-(8), (12), and (15): 

1.5/c (T e - Ti) = 1.5kT e {l + x e ) + (his - Xo ) + const - 

Setting {hv — %q ) ~ 1.5/cT e , and x^ ~ 1 — x e we have (T e — Ti) ~ T e . This means that, 
as the ionizing photons are absorbed in the stream flowing into the front and the gas is 
heated, the electron temperature is detached from the ion temperature. 

However, when the electron temperature exceeds 10 4 K, we can no longer neglect 
cooling of the gas by electron-collisional excitation of the L a line. Figure 1, which is based 
on our numerical calculations, shows that this effect limits the growth of T e , and that the 
energy exchange with heavy particles rapidly equalizes T e and Tj. It is important that 
the temperatures become virtually identical at a rather large optical depth: in the case 
shown in Fig. 1 (hv = 2 X q, 5 = 0.1, T = 5000 K, V = 300 kms" 1 ), T e differs from T { 
by less than 10% at r ~ 5. This implies that less than 1% of the energy pumped into 
the gas by the flux of ionizing radiation is radiated in the region where the flow has a 
substantially two-temperature character. Thus, this effect is small, and we will assume 
that the temperatures of the electrons and heavy particles upstream from the front are 
equal. The T(r) dependence for the model problem in the one-temperature approximation 
is also shown in Fig. 1. 

Based on the above arguments, we may use the relationships T e = Ti = T upstream 
from the shock front, and 

V \ k Tz [T(1 + Xe)] + VNkT{l + x ^ (jf) = Qph ~ QcolNe - (20) 
If we linearize equation (4) in accordance with (19), we obtain, denoting 3ft = k/ ^171^, 

sv=SN=SH ^mi±^M. (21) 

It is clear from this relationship that the value of SN does not exceed 1% for T < 2 x 10 4 K. 
In the same approximation, we can rewrite (20): 

Q 1 

-kV — [T(l + x e )\ = Q ph - Q col N e . (22) 

Heating and cooling of the gas upstream from the front are brought about by photons 
with energies hv\ > 13.6 eV (1 < / < N), which arrive from the postshock region. To 
reduce the required computer time, we divided the energy interval from 13.6 to 1000 eV 
into 50 intervals equal on a logarithmic scale, and treated each interval as a pseudoline 
with energy Ei equal to the mean value inside the interval and with flux equal to the total 
flux of the real lines and continuum falling in the given interval. We solved the above 
equations for these pseudolines, numbered in order of increasing energy. 
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We took the optical depth of the pseudoline with the maximum energy as the 
independent variable, in place of the z coordinate. Formally, this was done by a termwise 
division of all differential equations by equation (15b) with I = N. With this choice of 
independent variable, the equations describing the structure of the preshock region do 
not depend on iVo and H in the case of rarefied gas, and z(tn) linearly depends on N . 
At this stage, we did not consider the change of the preshock gas velocity due to gravity, 
and used a plane-parallel approximation, assuming that the extent of the preshock HII 
region was small (see below). All else aside, this makes it possible to reduce the number 
of free parameters. In the case considered of an intermediate-density gas, the structure 
of the preshock region should depend on N , since the level populations of some ions 
that appreciably contribute to the cooling function (e.g., C III, OIV, and O V) depend on 
density; however, this dependence proved to be rather weak. 

To solve the resulting system of equations, we must fix the value of V , which is a 
parameter, and specify the boundary conditions for T, z, 77, and x?. As stated above, the 
gas temperature far upstream from the front, i e., at z — >■ —00, should be of the order of 
several thousand Kelvin, and we treated T as another global parameter of the problem, 
varying it from 3000 to 6000 K. 

The boundary conditions for z and 77 (I < N) are obvious: z = T\ = for = 0. 
The question of the degree of ionization of the material far upstream from the front, 
i.e., specification of the boundary values for x°, is not so trivial. If we remain in the 
framework of a plane-parallel geometry, it is natural to assume that the gas is neutral at 
z — > —00. This automatically implies that the infalling matter is sufficiently abundant 
to completely absorb the ionizing radiation escaping from beneath the front. Generally 
speaking, however, this is not the case, since many CTTSs are X-ray sources. It is probable 
that, owing to the curvature of the magnetic lines of force, the motion of the infalling gas is 
appreciably non-rectilinear (radial), even at heights hx ~ R*/3 above the stellar surface. 
In this case, the hardest part of the postshock radiation will at some time cross the 
boundaries of the accretion column and will be able to reach the observer, in contrast to 
lower-energy photons, which are absorbed in the infalling gas at heights h < hx- (It is 
probably for this reason that the X-ray spectra of CTTSs resembles the radiation of a 
plasma with a temperature of ~ 10 7 K, while the temperature of the post-shock gas is 
considerably lower.) 

Thus, because of the curvature of the lines of force guiding the infall of material, some 
of the most energetic photons born at the point of intersection of the corresponding line 
of force with the stellar surface will not reach the gas at heights h > h X - However, this 
will to some extent be compensated by the absorption of hard radiation escaping from 
the limits of the accretion column above other points of the stellar surface. Since there are 
not currently any observational data about the field geometry in CTTSs, we will assume 
that, due to this compensation, the plane-parallel approximation remains valid at heights 
h > h x - More precisely, we assumed that when t n >> 1, H, He, N, O, and Ne are in atomic 
form, whereas C, Mg, Si, S, and Fe are singly ionized by Lyman lines of HI. Naturally, 
after solving the equations for this problem, we must verify that the extent of the preshock 
ionization zone does not exceed the radius R\ of the inner edge of the disk. 

In this approach, determining the heat and ionization balance of the preshock region 
is a boundary-value problem (some of the boundary conditions are specified at the shock 
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front and some at infinity), which we solved as follows. At sufficiently large values of 
t n — T~N ax , we set trial values of r™ a:E , (I < N — 1) and integrated the equations to 
Tn = 0. If, in this case, some T\ differed from zero by more than 0.01, we repeated the 
integration, changing the values of tJ 71 ^ accordingly. 

To solve this problem, we must have a set of values Si for I < N, which can be found 
only after determining the structure of the postshock region. However, for this, we must 
know the temperature and ion composition of the gas immediately upstream from the 
front. Thus, we obtained a self-consistent solution for the shock- wave structure using an 
iterative procedure: first, we determined the preshock gas parameters for a trial set of Si 
values, which made it possible to find first- approximation values of the Si; using these, 
we again performed the calculations for the preshock region, and so on, until the new Si 
values were virtually identical to the previous ones. 

4.2 The postshock region 

Let us now consider the simplifications we used to solve the complete system of equations 
for the post-shock region. First and foremost, we assumed that, at T > 10 4 K, the gas 
in the recombination region is transparent to the continuum and to its own "cooling" 
radiation; therefore, we did not take photoionization into account in this zone. In addition, 
we assumed that the density is so high (N e > 10 15 cm -3 in the transition region between 
the cooling zone and stellar photosphere, where the "cooling"-radiation photons are 
absorbed, that, due to thermalization, radiative cooling occurs mainly in the continuum 
at A > 912 A, to which all higher-lying layers are transparent. This assumption allows 
us to completely avoid considering the transition-region structure at this stage, since we 
are currently interested only in the line emission of "high-temperature" ions (see above). 
The structure and emission of the transition zone is important in analyses of the nature 
of the so-called "veiling" continuum of T Tauri stars and, perhaps, the hydrogen and He I 
emission-line spectra [30]. Note that the gas velocity in the transition region is nearly 
zero, whereas it is of the order of 300 kms -1 in the preshock region. Therefore, even if 
emission in Lyman-series lines carries away appreciable energy from the transition region, 
this should not strongly influence the thermal balance of the preshock zone. 

In the postshock region, it is more convenient to use the variable T s = T e x e + Tj in 
place of Tj. An equation for this variable can be derived by adding (6) and (7): 



Eliminating the derivative dN l /dz from this equation and equation (7), we obtain 
after some manipulation two equations, which we used in place of equations (6) and (7): 



Here, we have introduced the notation Cq = (T s ) 2 , C\ = 3tCo/V 2 2 , C 2 = 0.5V^ 2 /V 2 , v = 
V/V2, t s = T s /Cq, and y = T e x e /(T e x e ) 2 . As in (16), we denote parameters immediately 




V-2_ 2 dts_ ZflQcolXe 

N 2 V dz = f 2 kC 
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Vg 2 dy _ UeiXe - QcolXe {h + 2Ciy) 

N 2 V dz ~ 1.5JfeC 



(24) 
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downstream from the viscous jump with the subscript 2; fa = 3v 2 -2(1 + d + C 2 )v + Cit s , 
and / 2 = 3/i — 2Cit s . An additional relationship is provided by the algebraic equation 

v 3 - (1 + d + C 2 )t; 2 + Cit fl u + C 2 = 0, 

which can be obtained from (2), (4), and (5). 

It follows from (16) that, immediately downstream from the viscous jump, the ion 
temperature reaches T % 2 > 10 6 K, while the electron temperature only grows to ~ 3 x 
10 4 K. An energy exchange between the electrons and ions begins, and the temperatures 
are equalized at T™ ax ~ T 2 l /2. In this region, up to T e = 0.8T e maa: , we used T e x e as 
the independent variable; for this purpose, we performed a termwise division of all the 
differential equations by the equation for y. After this, we used T s as the independent 
variable, up until the time when the difference between T e and Tj became less than 1%. 
Finally, in the recombination zone, we used a one-temperature approximation and T e (l + 
x e ) as the independent variable, simultaneously leaving out equation (24) for T e . We solved 
the equations in this region up until the time when the gas temperature fell to T = 10 4 K. 

5 Discussion 

We solved the above system of radiative gas-dynamical equations numerically, using 
dedicated software for this purpose. We used a number of test problems to debug this 
program and verify its correctness, which ultimately gave us confidence in the reliability 
of the results obtained. In particular, the results of the numerical calculations coincide with 
analytical solutions we obtained for the case of a hydrogen plasma under certain additional 
assumptions. Our calculations of the ionization equilibrium for all elements included in 
the program coincide with the results of [36] in a (stationary) coronal approximation; 
however, the cooling function A(T) we obtained for this case differs from that of [36] on 
average by a factor of two, since we used more accurate formulas for the electron-collisional 
excitation rates, according to the recommendations of [37]. 

We used a special subroutine for the numerical solution of systems of stiff differential 
equations. To test its performance, we computed the ionization balance in a plasma 
with an arbitrary initial density of ions of various elements for the case of a stepwise 
change of its temperature: the solutions smoothly converged to the values obtained in the 
stationary case for a finite temperature. We also calculated the ionization balance and 
cooling function of an optically thin plasma for the case of isobaric cooling (this situation 
is similar to the behavior of the postshock gas [32]) and obtained agreement with [47], 
with appropriate variation of the atomic data. Finally, we note that, in the computations 
for the accretion shock wave, the sum of the x° values for each element differed from unity 
by less than 10~ 6 . 

Let us now turn to the results obtained. First and foremost, it follows from our 
calculations that the fluxes in the signal lines (see below) are virtually constant if 

(1) we take the initial temperature T of the infalling gas far upstream from the front 
to be in the interval from 3000 to 6000 K; 

(2) in our calculations of the preshock region structure, we do not take into account 
(pseudo) lines of the ionizing radiation, for which 5i < 10~ 4 ; 
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(3) we begin to integrate the equations describing the structure of the preshock region 
at values T^ ax exceeding 5 (in accordance with Section 2, we chose the optical depth of 
the line with the maximum energy of those with Si > 10~ 4 as the independent variable); 

(4) we vary the intensity of the z component of the chaotic magnetic field Hq frozen 
in the infalling plasma in the range from to 0.1 G. 

This last result deserves special attention. The chaotic magnetic field can suppress the 
electron thermal conductivity if the Larmor radius of an electron in this field is smaller 
than the mean free path of the electrons between collisions. Using formulas from [32], we 
can write this condition 



In the cases of interest for us, N < 10 13 cm -3 and T ~ 2 x 10 4 K in the preshock 
region (see below). Thus, the assertion in Section 4 a posteriori indicates the fundamental 
possibility of determining the structure of an accretion shock wave while simultaneously 
neglecting the electron thermal conductivity and the chaotic magnetic field frozen in the 
infalling plasma; this considerably simplifies the system of gas-dynamical equations. It is 
another question whether the accreted plasma is, in fact, magnetized, and whether we 
must in reality take into account the heating of the preshock gas by electron thermal 
conductivity; this question can be answered only after a comparison of the calculated and 
observed spectra. 

Below, we give a relatively detailed description of the shock-wave structure for the 
case when Vq = 300 kms -1 and iVo = 10 12 cm -3 , and then show how the main parameters 
of the wave change as functions of Vq and Nq. We will use 5000 K and 5 for the initial 
values of T and t^, respectively. 

Figure 2 shows the distribution of the electron and ton temperatures downstream and 
upstream from the shock front. We can see that the temperature of the gas approaching 
the front grows and reaches T\ = T(z — 0) ~ 1.96 x 10 4 K immediately upstream from 
the front. Since T\ < 2 x 10 4 K, in accordance with (21), we can consider our initial 
assumption that the variations of the preshock gas density and velocity are small to be 
justified. We will characterize the size of the preshock region by the distance from the 
front z pre at which the temperature of the infalling gas reaches a value exceeding T by 
50%; in the case considered, z pre ~ 1.8 x 10 9 cm. 

In accordance with (16), the ion temperature reaches ~ 2.66 x 10 6 K immediately 
downstream from the front, while the electron temperature only grows to T| ~ 4.9 x 10 4 K. 
Further, Tj monotonically decreases due to the exchange of energy with the electrons. 
The postshock electron-gas temperature first increases to T™ ax ~ 1.29 x 10 6 K at z T = 
z(T™ ax ) ~ 1.41 x 10 5 cm, then monotonically decreases as the radiative losses begin to 
dominate over the heating by the ions. At the peak of T e , the difference of the electron 
and ion temperatures is less than 1%, justifying use of a one-temperature approximation. 
The gas cools to T = 10 4 K at a distance z pst ~ 9.37 x 10 6 cm from the front. Here, we 
should note that both z pre <C R and z pst <C R*. 

Since Zt <C z pst and the width of the region in which the ions are heated is an order 
of magnitude smaller, our approximation of the front as an infinitesimally thin jump 
is justified. In addition, T™ ax differs from T™ ax /2 by only 3%; therefore, this suggests 
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Phc. 2: Variation of electron (T e ) and ion (T;) temperatures across the shock wave for 
V = 300 kms" 1 and N = 10 12 cm" 3 . 



that we can consider the entire postshock region in a one-temperature approximation 
(see the comment at the end of Section 3). However, there are two reasons why we did 
not do this. The first is associated with the signal-line profiles: in the one-temperature 
approximation, their widths immediately downstream from the front due to thermal 
motion will be systematically larger than observed. Second, in the region where the ion 
temperature exceeds T™ ax /2, nuclear reactions converting deuterium to 3 He can occur. 
We plan to treat the corresponding set of problems in a separate study, and here restrict 
our consideration to the comment that, according to our estimates, the energy release of 
these nuclear reactions is sufficiently small that they should exert virtually no effect on 
the shock-wave structure. 

Figure 3 shows the spectra of the ionizing radiation escaping from the gas cooling zone 
for three values of Vq and for logiVo = 11, or, more precisely, the values of Si for the 50 
energy intervals (pseudolines) . The dashed line at 5 = 10~ 4 enables us to determine which 
pseudolines play an appreciable role in the heat and ionization balance in the preshock 
region. In particular, it follows from this figure that, for Vq = 300 kms -1 , the line with 
the maximum energy is that with N = 44 (£44 ~ 571 eV), whose optical depth we used 
as the independent variable in this case. 

We will now consider how the gas ionization varies along the flow. The preshock H II 
region is approximately half the size of z we : zun — 7 x 10 8 cm. As the gas approaches the 
front, the degree of ionization of the hydrogen steadily increases, reaching = 7 x 10~ 4 . 
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Phc. 3: Emission spectrum generated behind the shock front for three values of the 
preshock gas velocity: V = 200 (a), 300 (b), and 400 kms -1 (c). In all cases, log N = 11. 



Immediately upstream from the front, 96% of the helium is in the form of He III; carbon 
and oxygen ions up to a charge of +4, nitrogen, neon, magnesium, and silicon ions up to a 
charge of +5, and sulfur and iron ions up to a charge of +6 all have abundances exceeding 
1%. 

Immediately downstream from the front, the time scale for heating of the electrons 
is much shorter than the time scale for ionization and recombination; therefore, the 
ionization equilibrium cannot be established after the rapid growth of the temperature. 
As a result, the abundance of, for example, CIV ions exceeds 1% up to T e ~ 750 000 K, 
and that of O VI ions becomes less than 1% only at T e > 1.1 x 10 6 K. Recall that, in an 
equilibrium situation, these ions are virtually absent at temperatures exceeding 2 x 10 5 
and 6 x 10 5 K, respectively. In short, the degree of ionization of the gas is considerably 
lower than the equilibrium value over nearly the entire region of growth of T e . 

As the ion and electron temperatures equalize, the time scale for the growth of T e 
decreases. Therefore, near the peak of the electron temperature, the ion abundances 
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roughly correspond to the equilibrium values for T e ma:E . However, soon after the electrons 
have begun to cool, the cooling time scale becomes much shorter than the recombination 
time scale, and the equilibrium again breaks down. Now the degree of ionization is higher 
than the equilibrium value: for instance, the CIV and OVI ions reach their maximum 
abundances only at temperatures of the order of 2 x 10 4 and 8 x 10 4 K, respectively! 

It follows from the above discussion, in particular, that the recombination of ions 
with charges less than six occurs at temperatures much lower than in an equilibrium 
situation. However, this means that, in our case, the role of dielectronic recombination 
is less important than usual. This suggests that, even if our description of dielectronic 
recombination is not entirely correct, this does not strongly affect the results obtained. 

We calculated the shock-wave structure in the approximation of volume radiative- 
energy losses; now we must test to what extent this is justified. We will start from the 
fact that typical values of the critical density for allowed dipole transitions are N cr = 
Aki/Cki ~ 10 17 — 10 18 cm~ 3 . For the part of the preshock HII region where N e ~ N , the 
optical depth in the L a line proved to be 10 6 (in the model with V = 300 kms -1 and 
iVo = 10 12 cm -3 ). Therefore, the probability of escape of a L a photon from this region, 
which is to order of magnitude equal to r" 1 / 2 ~ 10 -3 , considerably exceeds the probability 
of thermalization, equal to N e /N cr ~ 10 -6 . Beyond the region considered, the hydrogen 
becomes neutral, and the optical depth of the L a line grows rapidly, reaching ~ 10 8 at z — 
z pre ; however, the electron number density and gas temperature simultaneously decrease. 
This means that we can use the volume energy-loss approximation throughout the region 
where the contribution of the L a line to the cooling function plays an appreciable role. For 
resonant lines of other elements with A > 912 A, our assumption of the volume character 
of the preshock radiative losses is even better satisfied; for example, for the A1550 Aline 
of CIV, r ~ 6 x 10 3 . 

Similar estimates based on the results of our calculations show that the volume energy- 
loss approximation is also justified downstream from the shock front. In this very region, 
the gas turned out to be transparent to the Lyman continuum, confirming that we can 
neglect photoionization processes here. 

We have intentionally distinguished the postshock region with T > 10 4 K. The reason 
is that, at lower temperatures, HII and Hell, which should efficiently absorb photons with 
E > 13.6 eV traveling toward the star, begin to intensely recombine. In this transition 
zone between the shock wave and the stellar atmosphere, the gas velocity is nearly zero, 
and the gas density is several orders of magnitude higher than upstream from the front: for 
example, in the model considered, V ~ 0.9 kms" 1 and Af~4x 10 14 cm" 3 for T = 10 4 K. 
This indicates that we cannot use the volumeloss approximation for T > 10 4 K; instead, 
we must solve for the heating of the stellar atmosphere by external short- wavelength 
radiation taking into account the effects of radiative transfer in the lines and continuum, 
as was done in [48]. 

The calculations show that up to the time when our approach is no longer adequate, 
the abundances of ions with charges exceeding +1 become negligible for all the elements 
we have considered. It is likely that ions with charges +2 and higher will not appear in 
significant numbers in the transition zone, due to the high density of material there. This 
statement is based on the calculation of the transition-region structure carried out using 
the method we applied for the preshock zone, i.e., without taking into account the effects 



19 



Ta6jiim;a 1: Parameters of the structure of an accretion shock wave 



Vo, kms 1 


200 


300 


400 


log No, cm" 3 


12 


11 


12 


11 


Zpre; Cm 

znn, cm 

x(m) 

x(CW) 

pre 
T C IV 1548 

Ti, K 
T|, K 
T 2 \K 

rpmax TV" 
1 e i IY 

z(T ma x)i cm 
z pst , cm 


1.3 x 10 8 

1.3 x 10 8 
0.0027 
0.096 

2.2 x 10 2 
1.7 x 10 4 
4.2 x 10 4 
1.2 x 10 6 
5.6 x 10 5 
2.2 x 10 4 

1.4 x 10 6 


2.8 x 10 10 
2.1 x 10 10 
0.0008 
0.17 

5.6 x 10 3 
1.8 x 10 4 
4.5 x 10 4 

2.7 x 10 6 

1.3 x 10 6 

1.4 x 10 6 
9.3 x 10 7 


1.8 x 10 9 
7.0 x 10 8 

0.0007 
0.16 

5.9 x 10 3 
2.0 x 10 4 
4.9 x 10 4 
2.7 x 10 6 

1.3 x 10 6 

1.4 x 10 5 
9.4 x 10 6 


8.4 x 10 10 
2.4 x 10 10 
0.0005 
0.093 
2.0 x 10 4 
1.8 x 10 4 

4.6 x 10 4 

4.7 x 10 6 
2.3 x 10 6 
3.2 x 10 6 
2.0 x 10 8 



of thermalization of the radiation. The true temperature distribution inside the heating 
region may considerably differ from the calculated distribution, but we do not see why 
the efficiency of radiative recombination and charge exchange - which give rise to the 
near complete absence of highly- charged ions in our model of the transition zone - should 
decrease. 

Table 1 shows the main parameters characterizing the structure of an accretion shock 
wave as functions of Vq and N . We can see that the preshock gas temperature decreases 
somewhat with increasing N Q as a consequence of the reduced contribution of intercombination 
lines to the cooling function. This effect influences the geometric factors to a much lesser 
extent: z pre and z pst are nearly inversely proportional to N . The dependence of the 
parameters on Vo is rather trivial, and, in our view, does not require additional discussion. 

Table 2 lists the line fluxes for some ions, normalized to the total flux of the AA 1548, 
1551 doublet of CIV and expressed in percent. For determining the absolute flux values 
using relationship (18), the values of 5 [CIV 1550] are also indicated. Note that, in all our 
models, as in all T Tauri stars, the A1550 doublet of C IV proved to be the most intense 
line at wavelengths 1200-2000 A, apart from the L a A1215 A line, of course. 

Some ions have appreciable abundances in three regions: in the preshock HII region, 
immediately downstream from the front (where the electron temperature grows), and, 
finally, in the recombination zone. 

Therefore, generally speaking, these lines should have three-peaked profiles, due to 
the different gas velocities in these regions: Vo, Vb/4, and ~ 0. Note that multicomponent 
structure is, indeed, observed in some lines of highly-charged ions in the spectra of T Tauri 
stars [24]. Our calculations indicate that the relative intensities of these peaks change 
considerably depending on the values of Vo and A , and that this dependence is different 
for different ions. For instance, if log Ao = 12, the ratio of the fluxes of the CIV A1550 A 
doublet from the post-shock and preshock regions is 1 : 17 for Vo = 200 kms" 1 and 2 : 1 
for Vq = 300 kms" 1 . The flux from the post-shock region for the intercombination C III 
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Ta6jiim;a 2: Fluxes in lines of some ions (percent relative to the total flux in the CIV 
AA1548, 1553 doublet) 



Vo, kms 1 


200 


300 


400 


log N , cm~ 3 


12 


11 


12 


11 


o(CIV 1550) 


0.045 


0.042 


0.048 


0.049 


OVI 1032+1038 A 


900 


440 


330 


220 


NV 1239+1243 A 


54 


29 


22 


19 


Si IV 1394+1403 A 


16 


21 


26 


12 


0111] 1661+1666 A 


10 


92 


9.1 


69 


Si III] 1892 A 


2.2 


3.5 


0.9 


1.3 


CIII] 1908 A 


0.2 


1.6 


0.2 


0.8 



A1908 A line is negligible in all cases considered, while for the OVI A1031 A doublet, on 
the contrary, the intensity of emission from the preshock region is negligible. 

As we noted above, there is no thermalization of the signal-line emission. Since the 
parameters of the preshock H II regions obtained in our numerical calculations are nearly 
identical to those used in [30], we suggest that, if the preshock gas density does not 
strongly exceed 10 12 cm -3 , radiation at A > 912 A leaves the shock nearly unimpeded. 
This means that, in spite of the considerable optical depth, the intensity ratio for the 
components of the CIV A1550 doublet should be 2 : 1, since this is equal to the ratio of 
the electron-collisional excitation coefficients for the corresponding levels, which is nearly 
equal to the ratio of the statistical weights of these levels. The observed intensity ratio for 
the A1548 and A1551 lines in T Tauri stars is 2 : 1, but, in accordance with the arguments 
above, this should in no way be taken as proof that these lines are optically thin, as 
asserted in [3, 7, 24]. 

In the signal lines, whose optical depth considerably exceeds 1, the observed emission 
flux should depend on the area of the accretion zone and its orientation with respect to the 
observer, while the fluxes of optically thin lines are determined by the emission measure. 
Therefore, the relative intensities of the signal lines in the spectrum of a star, generally 
speaking, should differ from those listed in Table 2. For the same reason, in stars with 
different accretion-zone geometries, both the line profiles and the overall appearance of 
the emission-line spectrum should differ, even if their values of V and N are identical. 

6 Conclusion 

Based on the discussion in the previous section, we expect that, in the framework of 
our assumptions, we have correctly calculated the line-emission intensities for ions with 
charges above +1. We have also shown that the results of our numerical calculations 
confirm a posteriori the validity of our main assumptions, adopted to simplify the system 
of hydrodynamical equations used to model the structure of an accretion shock wave. It 
is important that the relative intensities and profiles of the spectral lines strongly depend 
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on Vo and Nq; consequently, the results obtained have diagnostic value. All this leads us 
to believe that comparison of the calculated fluxes and profiles of the C III, C IV, Si III, 
Si IV, etc., lines with those observed in the spectra of young stars may provide a base for 
elucidating the origin of the activity of these objects. 
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